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We study reaction-diffusion processes with concentration-dependent diffusivity. First, we de- 
termine the decay of the concentration in the single-species and two-species diffusion-controlled 
annihilation processes. We then consider two natural inhomogeneous realizations. The two-species 
annihilation process is investigated in the situation when the reactants are initially separated, namely 
each species occupies a half space. In particular, we determine the growth law of the width of the 
reaction zone. The single-species annihilation process is studied in the situation when the spatially 
localized source drives the system toward the non-equilibrium steady state. Finally we investigate 
a dissolution process with a localized source of diffusing atoms which react with initially present 
immobile atoms forming immobile molecules. 

PACS numbers: 05.40.-a, 82.20.-w, 66.10.C-, 05.70.Ln 



I. INTRODUCTION 

Random walks and their continuum descriptions in 
terms of convection-diffusion equations underlie numer- 
ous natural phenomena [TJ [5] . Our experience with these 
equations strongly influences our intuition. For instance, 
an instantaneous propagation of perturbations is a key 
property of the diffusion equation which appears puzzling 
when we learn it, but eventually we start to regard it as 
a general feature of parabolic partial differential equa- 
tions (PDEs) which distinguishes them from hyperbolic 
PDEs. It then comes as a surprise that the instantaneous 
propagation of perturbations does not generally apply to 
nonlinear parabolic PDEs, so that phenomenologically 
the predicted behaviors resemble the behaviors of hyper- 
bolic PDEs. Perhaps the simplest class of such nonlinear 
parabolic PDEs is 

^-V-( C Q Vc) (I) 

According to (JIJ , the spread of a localized initial profile 
into vacuum (that is, c(r,i — 0) = when r = |r| > R 
for some finite R) is not instantaneous whenever a > 0. 
More precisely, the front advances as r/ ~ ^i/(2+do) ( nere 
d is the spatial dimension); for r > r/, the medium is 
in the vacuum state [3H7]. Of course, there is no real 
contradiction with the standard lore as Eq. ([!]) is non- 
linear. Further, this spectacular effect occurs only when 
we consider the spread into vacuum (where the diffusion 
coefficient D = c a vanishes). 

Nonlinear parabolic PDEs form a fertile research area 
and the behavior of solutions to nonlinear parabolic 
PDEs is still poorly understood [7]. We emphasize 
that nonlinear parabolic PDEs similar to Eq. ([IJ and 
their microscopic brethren are by no means pathologi- 
cal. For instance, Navier-Stokes equations are nonlinear 
parabolic PDEs with temperature-dependent transport 
coefficients. Thus in describing heat conduction in the 
simplest case when the fluid velocity vanishes (heat con- 
duction without convection) we must solve a nonlinear 
parabolic PDE for the temperature T similar to Eq. ([I]) 



as the coefficient of thermal conductivity \ depends on 
the temperature [3l 0] . For instance, for the hard-sphere 
gas x indicating that for heat conduction in the 

hard-sphere gas is described by Eq. ([I]) with a = |. The 
concentration-dependent diffusion effect has been also 
observed in ionic crystals and oxides [SHllj. it is rele- 
vant in electrochemistry [321 [13] , and it has many other 
applications |14j . 

Nonlinear parabolic PDEs often arise as mathemati- 
cal models of reaction-diffusion processes |B] . Almost all 
studies of reaction-diffusion processes rely on the stan- 
dard diffusive transport. Reaction-diffusion processes 
with density-dependent hopping rates have appeared in 
a few studies, e.g., in Rcf. [15 annihilation processes 
have been investigated. Here we examine basic reaction- 
diffusion processes in the situation when the transport is 
described by (|lj) . 

In Sees. [Il}|lTl] we study the diffusion-controlled single- 
species and two-species annihilation processes, respec- 
tively. We examine the spatially homogeneous setting 
when particles are distributed at random with uniform 
concentration for the single-species annihilation process 
and with uniform and equal concentrations for the two- 
species annihilation process. In Sec. |IV| we consider the 
two-species annihilation process in the situation when 
the reactants are initially spatially separated, namely one 
half space is occupied by one species while the compli- 
mentary half space is occupied by another species. In 
Sec.[Vj we analyze the single-species annihilation process 
driven by a localized steady source. In Sec. |VI[ we study 
a dissolution process. This process involves three dis- 
tinct species, viz. diffusing atoms which is injected into 
a localized region and react with immobile atoms creat- 
ing immobile molecules. We summarized our results in 
Sec. En 

II. SINGLE-SPECIES ANNIHILATION 
PROCESS 

One of the simplest reaction-diffusion processes is the 
diffusion-controlled single-species annihilation. Symboli- 
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cally, this process is represented by the reaction scheme 



A + A 



(2) 



where identical diffusing particles (say atoms) are de- 
noted by A. The process ([2| postulates that when the 
two particles collide, they disappear. A true annihila- 
tion is occasionally possible (e.g., annihilation of domain 
walls in spin chains). In most situations, however, a col- 
lision between two atoms will lead to the formation of a 
diatomic molecule; if such molecules are stable (that is, 
they do not break back to atoms) and if molecules do 
not influence diffusing atoms, we can effectively ignore 
molecules and use the reaction scheme pi. 

The diffusion-controlled single-species annihilation 
process ([2| with concentration-independent diffusivity is 
described by the reaction-diffusion equation 



dc 
dt 



= DV 2 c - Kc 2 



(3) 



The reaction rate K depends on the details of the process; 
in the simplest case when spherical particles of equal radii 
R diffuse independently and immediately react upon col- 
liding, the reaction rate theory expresses K through the 
diffusion constant D and R, viz. 



K = IQnDR, 



(4) 



This well-known Smoluchowski formula [16-15] is valid 
in three dimensions. More generally, in d dimensions the 
reaction rate scales according to K ~ DR d ~ 2 [6]. In 
the spatially homogeneous setting, c = —Kc 2 , leading to 
the c ~ (Kt) -1 large time decay of the concentration. 
This is valid in d > 2 dimensions; for d < 2, the above 
formula predicts c ~ R 2 ~ d / (Dt), so the R— dependence 
is certainly wrong when d < 2. A heuristic way of un- 
derstanding the behavior in low dimensions relies on the 
basic properties of random walks and gives (in the spa- 
tially homogeneous setting) 



dc 
~dt 



(DRc 2 
[Dc 3 



Dc 2 [\n^\ 



d = 3 
d = 2 
d = l 



(5) 



from which we deduce well-known asymptotic behaviors 



t- L 

t- 1 \nt 

t -l/2 



d = 3 
d=2 
d= 1 



(6) 



In equation ([6| we have displayed only time depen- 
dence; using parameters D and R one easily restores the 
dimensionally-correct behavior, e.g., in two dimensions 
c - (Dt)- 1 \n(Dt/R 2 ). The asymptotic behaviors ^ 
have been established through a combination of simula- 
tions, heuristic arguments (as outlined above), and exact 
solutions in one dimension; they have been subsequently 
proven in Refs. [TU1 120] , 



The underlying microscopic process is simple to realize 
on a lattice. Each site is occupied by at most one parti- 
cle, and if a particle hops to the occupied site, both par- 
ticles annihilate. Particles undergo symmetric nearest- 
neighbor hoping. When the hoping rates are constant, 
the density decays according to Eq. (JgJ) . 

The simplest way to mimic the D = c° hoping rate is 
to postulate that the current hoping rate of each particle 
is inversely proportional to the distance £ between the 
particle and its nearest neighbor, namely 



hoping rate 



_ p-a/d 



(7) 



The quantity £~ x f d can be interpreted as a local con- 
centration, and hence the hopping rule ([7| mimics the 
macroscopic D ~ c a dependence. 

For the single-species annihilation process ^ with 
concentration-dependent hoping rate the generalization 
of ([5]) is very simple — one merely replaces the constant 
D by the concentration-dependent D — c a . This yields 



dc 
dt 



„2+a 



„2+a 
n 3+a 



d = 3 
d = 2 
d=l 



(8) 



A heuristic derivation of Q is equivalent to the deriva- 
tion of |5]). In three dimensions, one writes c = — Kc 2 , 
where the factor c 2 merely reflects the binary nature of 
the annihilation process; then the reaction rate theory 
(see [S] [T5J) gives K ~ DR oc c a thereby leading to 
c ~ — c 2+a as stated in In one and two dimensions, 
the recurrent nature of random walks becomes crucial. 
Generally we write c = — c/r, so we need to estimate the 
collision time r. In one dimension we use £ ~ %/ Dt where 
£ ~ c _1 is the typical separation between the nearest par- 
ticles. In two dimensions we need a slightly more compli- 
cated but still well-known formula Dr[ln ^j] ~ c~ l . 
The asymptotic behavior of the solutions to Eq. (JsJ) is 



ft- 1 /(i+«) d = 'S 

t -l/(l+a) pntjl/a+a) d = 2 
t -l/(2+a) d=1 



(9) 



For diffusion-controlled single-species annihilation pro- 
cesses with concentration-independent diffusivity there is 
a large body of knowledge (ranging from exact solutions 
in one dimension to proofs in all dimensions [T§1 I2TJ] ) 
supporting the asymptotic behaviors ^ and similar re- 
sults for more complicated processes including aggrega- 
tion, two-species annihilation, etc., see [2TM27] . There is 
no doubt in the validity of ([9| for single-species annihila- 
tion processes since the above (admittedly heuristic) ar- 
guments directly extend well-established results ([5])-([6]) 
to concentration-dependent hoping rates 

The generalization of rigorous work (T9l [20] to 
concentration-dependent hoping rates appears feasible, 
although it may require substantial effort. There is lit- 
tle hope of solving the model (even in one dimension in 
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FIG. 1: Particle positions in two-species annihilation in two 
dimensions. (The particle radii are enlarged for visibility). 



the simplest case of a = 1). Finally we note that for a 
ferromagnetic Ising spin chain supplemented with (con- 
servative) spin-exchange Kawasaki dynamics, the low- 
temperature behavior is well represented in terms of do- 
main walls undergoing single-species annihilation with 
diffusion coefficient proportional to their density [2"M3"U] . 
This corresponds to a = 1 and therefore according to 
Eq. ([9]) the density of domain walls should decay as i^ 1 / 3 
in one dimension. The celebrated i -1 / 3 decay law is very 
well-confirmed, both generally for conservative dynamics 
and specifically for the Kawasaki dynamics [28H30] , 



III. TWO-SPECIES ANNIHILATION PROCESS 

The two-species annihilation process is represented by 
the reaction scheme 



A + B 



(10) 



Here A and B denote diffusing particles of different type 
and the process (10) postulates that when the two dis- 



similar particles occupy the same lattice site, they imme- 
diately annihilate. The interesting situation arises when 
the initial densities are equal. (Otherwise, the minority 
species quickly disappears.) The critical dimension for 
this process is d c = 4. For d > d c , the mean-field de- 
scription is applicable, that is c ~ —DRc 2 ~ — c 2+a , and 
we recover the c ~ t _1 /( 1 + a ) decay. Below the critical 
dimension, the densities decay slower with time, namely 
as (Dt)~ d / 4 . This slow kinetics arises because opposite- 
species reactants organize into into coarsening domain 
mosaic (Fig. [I]) and annihilation can occur only along 
domain boundaries rather than throughout the system 
[23Tf27] . A heuristic explanation of the density decay is 
well-known and already appeared in textbooks [5J. This 
explanation is based on the domain picture. Namely, one 
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FIG. 2: Illustration of the three length scales in two-species 
annihilation in one dimension. 



argues that in a domain of size L ~ y/Dt one species was 
in majority (due to fluctuations in particle numbers), so 
the surplus survives and since it scales as \J coL d , the 
resulting density is c ~ L~ d \/ c§L d ~ *Jc$ (Dt)~ d ^. 

Replacing the hoping rate by its average value, D = c a , 
we obtain c ~ y/co (c a t)' d / i for d < 4. In the physically 
relevant dimensions, the asymptotic decay of the concen- 
tration is therefore 



c 2/(3a+4) t -l/(„+4/ 3 ) d = 3 

cJ /(o+2) t-Vta+2) d =2 
2/{a+4) t -l/ (a +4) d=l 



(11) 



Interestingly, in one and two spatial dimensions the 
coarsening domain mosaic (Fig. [T]) is characterized by 
three length scales [^. This is particularly obvious in 
one dimension (Fig. [2]) where one can identify the do- 
main size L ~ V Dt, the average spacing between adja- 
cent particles in the same domain £aa = @bb ~ c l , 
and the depletion zone between the domains. This latter 



quantity scales as £ab 



-1/4 



(Dt) 3 / 8 , see [31] or [6j. In 



the present case of the concentration-dependent hopping 
rate the three length scales characterizing the two-species 
annihilation process in one dimension are 



L ~ 

Iab 
Iaa 



2-a 



f2(a + 4) 



(12) 



(■BB 



Similarly generalizing the two-dimensional result |31j 
one gets 



L oc t"+ 2 



<-AB 



t-AA 



£ B B<xt'^+v (13) 



The validity of these results is again very plausible as they 
directly generalize well-known results for the standard 
diffusive transport. Numerical verifications could be very 
challenging even in one dimension. Indeed, if the hopping 
rate is proportional to the concentration (a = 1), the 
length scales (12) become L oc i 2 / 5 , Iab oc t 3 / 10 and 
Iaa t 1 / 5 , so the exponents differ only by 0.1 which is 
hard to accurately measure. 



IV. INHOMOGENEOUS TWO-SPECIES 
ANNIHILATION 

In this section we return to the two-species annihilation 
process from Sec. |III[ but instead of uniform initial den- 
sities we assume that the reactants are initially spatially 



4 



/ 


1 


— — / 

\ 1 
\ 1 

\^ / 


\ 


-* / j 


c > x 



w 



FIG. 3: Sketch of the concentration profiles (solid curves) and 
the reaction rate Kcacb (dashed, and considerably magni- 
fied) in two-species annihilation with initially separated com- 
ponents. The depletion zone width is I ~ y/t, while the reac- 
tion zone width is w ~ t 1/{4a+6) . 



separated. Specifically, we assume that the right (resp. 
left) half space is initially occupied by A (resp. B) par- 
ticles (Fig. [3]) . Physically it may be realized by having a 
membrane at x — separating the reactants, removing 
the membrane at time t = 0, and observing the subse- 
quent reaction kinetics. The governing reaction-diffusion 
equations read 







Dx V dx ) 



dcA 
dt 

dc B _ d ( d 9c b\ 
dt dx\ dx J 



- Kcacb 

- Kcacb 



(14a) 
(14b) 



where ca = ca(x, f) and cb = Cb{x, t) denote the concen- 
tration of each species at position x at time t. We assume 
that diffusion coefficients of both species are equal and 
depend on the total concentration ca + Cb. 
The initial densities are 



c A (x,t = 0) 



Co, x > 0; 
0, x<0; 



and cs(x,t — 0) = ca{— x,t = 0). In the constant- 
diffusion case, this problem has been originally investi- 
gated in Refs. [321133]; various generalizations have been 
treated later, see [3~4Tf4"0] and references therein. 

We now show how to handle the concentration- 
dependent case; in particular, we determine the width 



w of the reaction zone. Subtracting (14b) from (14a) we 
find that c : 



ca 



Cb satisfies a diffusion equation 




dc 
dt 



dx\ dxJ 



(15) 



The domain of one species acts as a nearly fixed absorb- 
ing boundary condition for the opposite species. Indeed, 
the ratio w jl of the width of the reaction zone w to the 
width of the depletion zone I asymptotically approaches 
to zero as we shall see below, and therefore in analysis 
of Eq. ( |15[ ) we can treat the reaction zone as an inter- 
face and summarize its influence by using the absorbing 
boundary condition. Thus the density profile of each par- 
ticle species, namely A particles in the half space x > 



and B particles in the half space x < 0, is close to that 
of independent diffusing particles in the presence of an 
absorbing boundary at x — 0. The diffusion coefficient 



depends on both concentrations, D = (ca + cbY 1 , so (15) 
is not a closed equation. Away from the reaction zone, 
however, one of the species dominates. Thus e.g. for 
x 3> w we can re- write ( 15 1 as 



c t = (c a c x ) 2 



(16) 



We want to solve ( 16 ) in the half-space x > subject to 
the initial condition 



c(x > 0,i = 0) = c 
and the absorbing boundary condition 
c(x = 0,t > 0) = 



(17) 



(18) 



The problem ( 16 )-( 18 1 admits a self-similar solution 
c{x,t) = coF({), £ = ^L= (19) 



'r a t 



The scaling function F(£) is a solution of the differential 
equation 



(F a F')' + §£F' = 
subject to the boundary conditions 

F(0) = 0, F(oo) = 1 



(20) 



(21) 



The boundary- value problem ( 20 )— ( 21 ) does not admit 
an analytical solution, but we are interested in the small 
£ behavior which can be extracted analytically (up to 
numerical factors). Indeed, an asymptotic analysis gives 



F~£ 



l/(a+l) 



when £ — » 



(22) 



We now estimate the width w of the reaction zone by 
balancing the diffusive flux into this zone with the total 
rate at which particles are annihilated. The flux is given 
by 2|c a C;r| and using Eqs. (19) and (22) we find that that 
the flux scales as c^ +1 / 



Using the reaction-diffusion 
equations (|14a|)-(|14b|, the number of reactions per unit 



cgt. 



time equals the integral of Kcacb over the extent of the 
reaction zone. We estimate this integral as c a+2 w (since 
K ~ c a ) and use c = c(x = w,t) as the typical concen- 
tration. We simplify c a+2 w using Eqs. (19) and p2| and 
equate the result with the flux to yield 



ill 

.,+1 



from which 



wc. 



w 



a+2 



/Cgt 



/cgt 



_-(o+2)/(4o+6) ,l/(4a+6) 
c C 



(23) 



We see that the width of the reaction reaction zone ex- 
hibits a rather slow growth w ~ /; 1 /( 4a + 6 ) ; while the width 
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of the depletion zone grows diffusively, I ~ \ft. Hence the 
ratio w/£ ~ £-(a+i)/( 2o + 3 ) asymptotically vanishes and 
this provides an a posteriori justification of the usage of 
the absorbing boundary condition (18 1. 



The above analysis applies when the annihilation pro- 
cess is happening in three (or higher) dimensions. We 
now show how to handle one- and two-dimensional cases. 
(See Rcfs. [35 37J for the analysis of such low-dimensional 
settings in the case of constant diffusion coefficient.) In 
these situations, we can use ( 16 )— ( 22 1 , so the flux still 



scales as Cq 



I cgi. The reaction rate c a+2 should be re- 



placed by c a+3 in one dimension, and by c a+2 [ 
in two dimensions, see Eq. Using these results we 
estimate the total rate at which particles are annihilated 
and balance the corresponding estimates with the flux. 
This yields 



W 



-1/2 ,l/(2a+4) 
L 1 



(24) 



in one dimension and 

W ~ c H Q + 2 )/( 4a + 6 ) £l/(4a+6) (lnt) (a+l)/(2a+3) (35) 



in two dimensions. Combining (|23|)— ( 25 ) and ignoring 



the dependence on the initial concentration we conclude 
that the width of the reaction reaction zone grows as 



t l/(4a+6) d = s 

f l/(4a+6) Q n ^(a+l)/(2a+3) d _ 2 
t l/(2a+4) d = 1 



(26) 



V. 



ANNIHILATION PROCESS DRIVEN BY A 
LOCALIZED SOURCE 



We now return to the single-species annihilation pro- 
cess and consider the system which is initially empty 
and is driven by a localized source, say at the origin, 
which is turned on at time t = 0. In the case of classi- 
cal (concentration-independent) diffusion, the qualitative 
behaviors are largely understood [411 142] ; we now review 
the approach and then apply it to non-linear diffusion. 



Constant Diffusion 



The corresponding reaction-diffusion equation reads 



dc 
Of 



= DV 2 c- Kc 2 + J5(r)9(t) 



(27) 



Here D is the constant diffusion coefficient, J is the 
strength of the source, and 9{t) is the Heaviside step 
function assuring that the source is turned on at t = 0. 
Equation (27) is applicable when d > 2, while for d = 2 
and d = 1 the reaction term must be modified according 
to the right-hand side of ([5]). 



Assuming the emergence of the stationary concentra- 
tion profile one gets 

= DV 2 c - Kc 2 + J 8(r) (28) 

Balancing the first two terms, DV 2 c = Kc 2 , leads to 



(4-d) 



2D 1 

~K r 2 



(29) 



Thus we see the emergence of the upper critical dimen- 
sion d c — 4. Indeed, the result of Eq. ( 29 ) clearly holds 
only when d < d c = 4. For d > d c = 4 the reac- 
tion term becomes irrelevant far away from the source. 
Therefore one must merely solve the Laplace equation 
DV 2 c + J 5(r) = 0. The solution is 



J 1 

£) j.d—2 



(30) 



At d = d c it is natural to expect logarithmic corrections. 
Indeed, one finds El SH [42] 



2D 

~K r 2 



1 



lnr 



(31) 



The behavior ( 29 1 is actually valid as long as the dimen- 
sion exceeds the critical d c = 2. At d = d c we must use a 
modified reaction term, so we balance V 2 c ~ c 2 /ln(l/c) 
to yield c ~ r~ 2 lnr. Finally for d = 1 we balance 



and get c 



(29H31) we arrive at [BHD 12] 



Combing these results with 



c(r) 



r -{.d-2) d > 4 
r'^lnr)- 1 d = 4 
r" 2 d = 3 

d = 2 
d=l 



(32) 



r lnr 



It is interesting to estimate the asymptotic growth of 
the total number of particles in the system. In principle, 



N{t) 



drQ, d r d - l c{r,t) 



where fi^ is the surface area of unit sphere in the d dimen- 
sional space: Oi = 2, VI2 — 2ir, = 4ir, etc. The con- 
centration is close to the steady-state distribution (32) 
for r < R(t), and it quickly vanishes when r > JUT). 
Therefore 



N(t) 



K(t) 



drQ d r d - l c{r) 



(33) 



In the case of constant diffusivity the radius of the occu- 
pied sphere grows as 



(34) 



31(f) - VDt 



G 



Inserting d32l) and pS) into Q one finds SB H2] 



Finally plugging (37) into (33 1 and using (38) we obtain 



N(t) 



For d > 4, there is so much room that particles do not 
'see' each other once they are sufficiently far from the 
origin. Overall the finite fraction of particles survive; the 
total number of surviving particles is of course smaller 
than the total number of introduced particles Jt, but it 
still grows linearly. For d < 4, the vanishing fraction of 
the introduced particles survives. 



't 


d > 4 






't 


d> d c 


t(lnt)- 1 


d = 4 








d = d c 


< t 1 ' 2 


d = 3 


(35) 


N(t) ~ < 


t (d-2) / (2+2a) 


2<d<d c 


(hit) 2 


d = 2 






(lnt) 2 


d = 2 


\nt 


d = l 






Int 


d = l 



In particular, in the most relevant case of d 
number of particles grows according to 

N(t) ~ tV(2+2a) 



(39) 



3 the total 



(40) 



Note also that for d = 1,2 the growth of the total num- 
ber of particles is very slow, namely logarithmic; interest- 
ingly, the asymptotic growth laws for the total number of 
particles are independent on a in one and two dimensions. 



B. Concentration-Dependent Diffusion 

In the case of concentration-dependent diffusion, the 
governing equation reads 



dc 
di 



= V- (c a Vc)-c 2+a + JS(r)9(t) 



(36) 



when d > 2, while for d = 2 and d — 1 the reaction term 
must be modified according to the right-hand side of pi . 
We again assume the emergence of the stationary con- 
centration profile. Balancing the reaction and diffusion 
terms V • (c Q Vc) = c 2+a one finds c = 2(4 + 2a - d)r~ 2 
which tells us that the upper critical dimension is now 
given by d c = 4 + 2a and the above expression for the 
steady-state concentration is valid when 2 < d < d c . Re- 
peating the same arguments as above we determine the 
concentration in every dimension to yield 



c(r) 



■ r -(d-2)/{l+a) d>d c =4: 

-2n_w d = d c 

2<d<d c 



2a 



r "(lnr) 

2 



r 

r~ 2 lnr 



(37) 



d= 2 
d=l 



C. Examples 

The simplest example of the concentration-dependent 
diffusion is D = c. This is relevant e.g. for the heat 
transfer in the so-called Maxwell gas [43 45 since the co- 
efficient of thermal conductivity (and generally transport 
coefficients in the Maxwell gas) is proportional to tem- 
perature. In this case (38) shows that in the physically 
interesting dimensions d = 1, 2, 3 the front propagates as 



<V4 d = 3 

W) - { t 1/4 (lni) 1/4 d = 2 
t 1 / 3 d = 1 



(41) 



Another interesting example corresponds to D = s/c. 
(For the hard-sphere gas, the transport coefficients are 
proportional to the square root of temperature.) In this 
case the front propagates as 



<V3 d = 3 

&(t) - { iVSQntji/s d = 2 
i 2 / 5 d = 1 



(42) 



In the classical case, the position of the front (34) 
was universal (independent on the spatial dimensional- 
ity). Now the situation is different. To determine the 
time dependence of the front position we use the chief 
formula 51 - VDt together with D - [c(3?)] a . Thus we 
have Oi 2 ~ [c(3i)] a t, and using (37) to estimate c(3?) we 
conclude that 



D. 



( t (l+a)/(2+da) d>d c 

t l/(2+2a)( lnt ya/(2+2a) d = 

t V(2+2a) 2<d<d c 

t l/(2+2a)( m ^a/(2+2a) rf=2 

*l/(2+a) d = ! 



(38) 



Exact Solution of the Mean-Field Equations in 
One Dimension 



are mean- 
?ly provide 



Reaction-diffusion equations ([3]), ([27]), (36 
field in nature. These equations presuma 
asymptotically exact descriptions in three dimensions. In 
two dimensions, the mean-field treatment of the annihi- 
lation process is slightly wrong (the discrepancy is loga- 
rithmic), and one dimension the mean-field approach is 
flawed. Even when reaction terms are chosen to assure 
the correct qualitative behaviors, there is no closed equa- 
tion for the density when the spatial dimension is below 
critical. In one dimension an exact description is feasible 
only in the case of constant diffusion coefficient [4"T1 146] . 
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It seems impossible to generalize these results to the case 
of the concentration-dependent diffusion. 

Nevertheless, it is methodologically interesting to solve 
the mean-field equations in the one-dimensional setting. 
Thus we want to solve an equation 



dc 
dt 



dx 



dc 
dx 



„3+o 



J6(x)6(t) (43) 



This is a rather complicated initial-value problem. We 
limit ourselves with the steady-state solution. The gov- 
erning ordinary differential equation 



d 

dx 



dx J 



c 3+a + J5(x) = 



(44) 



is solvable. Equation (c a c')' 



„3+a 



where prime de- 



notes the derivative with respect to x, admits an exact 
solution c = A/x with A — \fd? — y/4 + 2a. The transla- 
tional invariance of equation (c a c')' = c 3+a implies that 
its general solution is given by c = Aj(x + Xq), where Xq 



is an arbitrary constant. Returning to (44) and invoking 
the x ■<-> — x symmetry we see that the solution must be 



.4 



x 



x 



2A 1 



J 



(45) 



The constant xq in ( 45 ) was fixed by integrating ( 44 ) over 



a small region around the origin which gives 



dc 
dx 



x=+0 



x=-0 



-J 



and allows to express xq via the strength of the source. 

Solutions similar to ( 45 ) have been obtained in [B] [5T] . 
A general case with reaction term c n has been recently 
investigated in Ref. |15j where it was additionally demon- 
strated that some of these solutions fit experimental 
data [47] for morphogen gradient formation. (Morphogcn 
gradients play a crucial role in developmental biology 
[4*51 14*9"] . e.g., they appear to be precursors to cell dif- 
ferentiation [50] 151] .) 



VI. DISSOLUTION PROCESS 

Here we consider another reaction-diffusion process, a 
dissolution process, which involves three distinct species. 
One species is composed of diffusing particles, while the 
particles constituting two other species are immobile. 
Diffusing atoms (species A) are injected into a small lo- 
calized region of a d— dimensional lattice. The entire lat- 
tice is initially occupied by immobile atoms (species B), 
one B atom per lattice site. Whenever an A atom hops 
to a lattice site occupied by a B atom, two atoms re- 
act to form an inert stable molecule (species B*). This 
dissolution process is described by the reaction scheme 



A(diffusing) + S(substrate) B* (stable) 



(46) 



This model mimics a number of important industrial 
chemical processes including dissolution of solids [52] . 
electropolishing [53], corrosion and etching [54TI56] . and 
erosion [57] . Here we consider the dissolution process 
with concentration-dependent diffusion coefficient. 

The reaction ( 46 ) proceeds at a certain finite rate. In 



many applications this rate greatly exceeds the hopping 
rate. Thus we shall assume that the reaction proceeds 
instantaneously, so each site contains either a B atom or a 
B* molecule. As the process develops, the system can be 
separated onto two parts: The droplet that contains no 
B atoms (every lattice site inside the droplet is occupied 
by a B* molecule and can contain an arbitrary number of 
A atoms) and the rest of the system that contains only B 
atoms. The droplet is asymptotically a ball of radius R(t) 
and its boundary moves more and more deterministically 
when t> 1. We shall ignore fluctuations, so the results 
should be applicable for sufficiently large time. 

The concentration c(r, t) of A atoms (that is, the av- 
erage number of A atoms per lattice site) satisfies a non- 
linear diffusion equation 



dc 
di 



= D 



dr 2 



d-1 d 
r dr 



„l+o 



1 



J6(r)9(t) (47) 



inside the droplet < r < R(t). The diffusion equation 
( 47 ) should be supplemented by the adsorbing condition 



c(r = R(t),t> 0) = 
and the Stefan boundary condition 



d_R = Dca dc 
dt dr 



r=R 



(48) 



(49) 



The boundary moves and its position has to be deter- 
mined in the process of solution. Therefore mathemati- 
cally we arrive at the Stefan problem 58 , 59J . We analyze 
( 47 )-( 49 ) using the same approach [55J [BD] as in the case 
of the concentration-independent diffusion coefficient. 

The original process occurs on the lattice and there- 
fore we set the lattice spacing to unity; this implies that 
the spatial coordinates r, the droplet radius R(t), and 
the concentration c(r, t) are all dimensionless quantities, 
while the hopping rate D and the strength of the source 
have the dimension of inverse time: [D] = [J] = l/(time). 
The ratio J/D, the dimensionless flux, plays an impor- 
tant role in the problem. 

Consider for concreteness the most important in ap- 
plications [5~2] [53] [57] two-dimensional setting. Seeking 
solution in the scaling form 



c(r,t) = c(0, £ 



r 
R 



we reduce the diffusion equation (471 to 



A2 



RR £c' 
D c a 



(50) 



(51) 
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where (•)' = d(-)/d£ and R = dR/dt. Equation (51) is 



consistent if RR/D, which is in principle a function of 
time, is actually a constant. Hence we write 



R 2 = 2(3Dt, 



and recast (51 1 into 



c " + r y + a (c ')2 + ^ 

c c a 







which must be solved subject to 

c = 0, c a c + (3 = when £ = 1 



(52) 



(53) 



(54) 



Equation (54} implies that c~ \J3(1 + a)(l - 0] 1/(1+a) in 
the £ — ?> 1 limit. The conservation of the total number of 
A atoms gives 



Jt = I c(r, t) 2?rr dr + 7T.R 2 
Jo 



(55) 



Using Eqs. (50) and (52), we re- write (55) in the form 
J 



(56) 



which implicitly determines /3 in terms of the dimension- 
less flux J/D. The right-hand side of (56) depends of 
course on /3 since the concentration c(£) is determined 
by solving Eqs. (53)-(54) which contain [3. 



VII. SUMMARY 

In annihilation processes the concentration decays to 
zero, and even in annihilation processes driven by a lo- 
calized source the concentration asymptotically vanishes 
when the distance from the source increases to infinity. 
Hence it is vitally important to know the behavior of D(c) 
in the c — > limit. If D(0) = 0, the behaviors qualita- 
tively differ from the classical case of constant diffusivity. 

We studied single-species and two-species diffusion- 
controlled annihilation processes in the situation when 



the diffusion coefficient vanishes algebraically, D ~ c a , in 
the c — > limit. In the homogeneous case, the critical di- 
mensions are d c = 2 (for the single-species annihilation) 
and d c = 4 (for the two-species annihilation), and in the 
physically interesting dimensions the long-time asymp- 
totics are respectively given by ([£]) and |Tl| ). For the 
two-species annihilation we also investigated the width of 
the reaction zones, both in the homogeneous setting (the 
width of the zone between adjacent domains) and the 
inhomogeneous setting (when initially the reactants oc- 
cupy complimentary half spaces); chief results are given 
by Eqs. ([l2l-(fl3]) and (p6|. We then studied the single- 



species annihilation process in the situation when the 
spatially localized source drives the system toward the 
non-equilibrium steady state. We showed that two crit- 
ical dimensions demarcate different behaviors: d c = 2 
which coincides with the critical dimension of the homo- 
geneous process, and d c = 4 + 2a which depends on the 
exponent a. Our main results are given by (37)-(39). 



We also investigated the dissolution process involv- 
ing diffusing and two immobile species of particles. As 
the process develops, the system can be separated onto 
the spherical droplet that contains no original immo- 
bile atoms and the untransformed rest of the system. 
Mathematically, one needs to solve the Stefan problem 
as the radius of the droplet is determined in the pro- 
cess of solution. In the most important in applications 
two-dimensional setting, the solution is self-similar and 
this allowed us to reduce the problem to an ordinary dif- 
ferential equation. This equation is analytically soluble 
only in the case of concentration-independent diffusiv- 
ity when the governing equation is linear [55]. In this 
classical case, the probability that a diffusing atom has 
not become a part of a molecule has been recently de- 
termined [60] . It appears possible to generalize analysis 
of Rcf. [60] and probe first passage characteristics in the 
case of concentration-dependent diffusivity. 

Finally we mention recent studies [5TJ [S5] of fluctua- 
tions in diffusion processes with concentration-dependent 
diffusivity. An interesting challenge is to probe the 
role of fluctuations for annihilation processes with 
concentration-dependent diffusivity. 
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